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Abstract 

Hadronic annihilation decay rates of P-wave heavy quarkonia are given to next- 
to-leading order in both ag and v‘^. They include ten nonperturbative parameters, 
which can be rigorously defined as the matrix elements of color-singlet and color- 
octet operators in NRQCD. We expect these papameters will be determined from 
lattice calculations in future. 

PACS number(s): 13.25.Gv, 12.38.Bx 

BBL (Bodwin, Braaten and Lepage) |l[] factorization formalism has provided a rigorons QCD 
prediction for heavy qnarkonia decay and prodnction. In this formalism, the decay rate can be 
written in terms of two expansions in strong conpling constant and the relative velocity of 
qnark and antiqnark in qnarkoninm. Nonrelativistic results including QCD radiactive corrections 
have been given for S-wave and P-wave decay rates. However, in heavy qnarkoninm system, 
especially for charmonium, relativistic effects are very important. In cc system, ~ 0.2, and is 
comparable with as{mc){^ 0.3), therefore relativistic corrections can not be neglected. The case 
of S-wave have been studied in |^. In this paper we will consider the relativistic corrections for 
P-wave decay rates and give complete expressions, which are accurate to next-to-leading order 
in both as and 

In NRQCD, the effect of annihilation can be taken into account by adding 4-fermion oper¬ 
ators to NRQCD Lagrangian: 
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where the sum is over all possible local 4-fermion operators On that annihilate and create a QQ 
pair, and is the scaling dimension of On- The short distance coefficients fn{o(s) can be com¬ 
puted by matching perturbative amplitudes for QQ scattering in NRQCD with the corresponding 
amplitudes in full QCD. The annihilation rate of a quarkonium state H to light hadrons (LH) 
can be written as 


T{H —> LH) = 2Im < H\5Ci_ferm.ion\H > (2) 

At a given order in the number of matrix elements can be reduced to a finite number by 
using velocity scaling rules for the matrix elements . These scaling rules consist of that for the 
operators and for the probabilities of the Fock states that give the leading contributions to the 
matrix elements. 

To next-to-leading order in the decay width of four P-wave quarkonium states can be 
written as 
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Oi{^Po) = D) ■ D) ■ atlj 

Oi{^Pi) = D xa)x ■ X+(-^ D xa)^ 

Oi(3p2) = ^+(-1 

OiCPi) = ^+(-| D)x ■ x+(-^ D)V^ 

OsC^o) = 

08('^i)=V^+TVx-X+TV^ 

^i('^o) = D) ■ axx^{~ D) ■ + /i.c] 

B) VVx+(-^ B)(V^’V + h.c] 

V.QSo) = + h.c] 

Vs{^S,) = 1[V;+TV(-^ ^fx ■ X+TV^ + /i.c] (7) 


where D is the space component of covariant derivate %[) and x are two component operators 
of quark and antiquark respectively. Coefficients are 
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where 

i >0 = - 2n/), 

andGF = ^, C^ = A'o. 

Comparing with the leading order results, we have added four new hrst order operators 
Pi(^Fo), ’Pi(^-P 2 ), Vsi^So), and Vsi^Si), which coefficients are calculated only to leading order 
in ag. The coefficients of zeroth order operators have been given in IQ. Therefore we neglect 
higher order terms such as a^T, v^T and a^u^T, and only keep those accuracy to hrst order in 
as or Since the hrst order relativistic corrections have been involved, there are no dehnite 
relations for matrix elements of zeroth order operators Oi{^Pj), Cl8(^'S’i) (J=l,2,3), Cli(^Pi) and 
08 (^ 50 ). But for hrst order operators, due to heavy quark spin symmetry, we have 

< Xol'Pi(’-Po)IXi) > = < X2|'Pl(’-P2)|X2 >, (18) 

< > = < h\V,{S)\h > . (19) 

These two matrix elements together with eight zeroth order operators’ matrix elements consist 
of ten nonperturbative parameters. We know that there are only two parameters Hi and Hg 
in the case of nonrelativistic limit, which can be phenomenologically determined from two 
experimental data. Now we have not enough experimental values to determine all these ten 
parameters. But their rigorous dehnition has been given in NRQCD, and we expect lattice 
calculations can provide their numerical results. 

Now we give a reasonable estimate for the theoretical errors in our expression of decay rates. 
The three main sources of theoretical error are the neglecting terms a^T, u^T, and agv'^T, which 
can be estimated to be 9%, 4% and 6% respectively for charmonium system. Combinig the three 
errors, we obtain that the theoretical error is about 12% by using the standard formulas for 
propagating independent errors. For bottonium system, the error is much small, therefore our 
formulas (1), (^), and @ can give strong theoretical predictions for P-wave quarkonia decay. 
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